A PROOF OF THE SMOOTHNESS OF THE FINITE TIME HORIZON AMERICAN 

PUT OPTION FOR JUMP DIFFUSIONS 



ERHAN BAYRAKTAR 

Abstract. We give a new proof of the fact that the value function of the finite time horizon American put 
option for a jump diffusion, when the jumps are from a compound Poisson process, is the classical solution 
of a free boundary equation. We also show that the value function is C 1 across the optimal stopping 
boundary. Our proof, which only uses the classical theory of parabolic partial differential equations of 
[7, 8], is an alternative to the proof that uses the theory of vicosity solutions (see [14]). This new proof 
relies on constructing a monotonous sequence of functions, each of which is a value function of an optimal 
stopping problem for a geometric Brownian motion, converging to the value function of the American put 
option for the jump diffusion uniformly and exponentially fast. This sequence is constructed by iterating 
a functional operator that maps a certain class of convex functions to classical solutions of corresponding 
free boundary equations. On the other hand, since the approximating sequence converges to the value 
function exponentially fast, it naturally leads to a good numerical scheme. 



1. Introduction 

Let (f2, J 7 , P) be a complete probability space hosting a Wiener process W = {Wt,t > 0} and a Poisson 
random measure N on K + x IR + , with mean measure \v{dx)dt (in which v is a probability measure on 
independent of the Wiener process. We will consider a Markov process S = {St',t > 0} of the form 



dS t = yiS t dt + aS t dWt + S t ^ [ (z - l)N(dt, dz). (1.1) 



In this model, if the stock price jumps at time t, then it moves from St- to St = ZSt~, in which Z is a 
positive random variable whose distribution is given by v. Note that when Z < 1 the stock price jumps 
down and when Z > 1 the stock price jumps up. In the Merton jump diffusion model Z = exp(Y), in 
which Y is a Gaussian random variable. We will take /i = r + A — A£, in which ^ = xv(dx) < oo, so 

that (e~ rt St)t>o is a martingale, i.e., P is a risk neutral measure. The constant r > is the interest rate, 
and the constant a > is the volatility. We assume the risk neutral pricing measure P, and hence the 
parameters of the problem, are fixed as a result of a calibration to historical data. The value function of 
the American put option pricing problem is 

V(x,T):= sup E x {e- TT {K - S T ) + }, (1.2) 

tG<S ,t 

in which 5o,t is the set of stopping times (of the filtration generated by W and N) that take values in 
[0,T], and E x is the expectation under the probability measure P, given that Sq = x. 
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We will show that V is the classical solution of a free boundary equation and that it satisfies the 
smooth fit principle, i.e., V is continuously differentiable with respect to its first variable at the optimal 
stopping boundary. We argue these facts by showing that V is the fixed point of an operator, which we 
will denote by J, that maps a given function to the value function of an optimal stopping problem for 
a geometric Brownian motion. This operator acts as a regularizer: As soon as the given function / has 
some certain regularity properties, we show that J f is the unique classical solution of a corresponding free 
boundary equation and that it satisfies the smooth fit principle. The proof of the main result concludes 
once we show that V has these certain regularity properties. In this last step we make use of a sequence 
(which is constructed by iterating J starting with the pay-off function of the put option) that converges 
to V uniformly and exponentially fast. Incidentally, this sequence yields a numerical procedure, whose 
accuracy versus speed characteristics can be controlled. Each element of this sequence is an optimal 
stopping problems for geometric Brownian motion and can be readily calculated using classical finite 
difference methods (see e.g. [18] for the implementation of these methods). An alternative proof of the 
regularity of V was given in [14]. This proof used a combination of the results in [8] and the theory 
of viscosity solutions. In particular the proof of Proposition 3.1 in [14] is carried out (details are not 
provided but hinted) using arguments similar to those used in the proof of Proposition 5.3 in [15]. The 
latter proof uses the uniqueness results of [9] for viscosity solutions. 

The infinite horizon American put option for jump diffusions were analyzed in [3] using the iterative 
scheme we describe here. The main technical difficulty in the current paper stems from the fact that 
each element in the approximating sequence solves a parabolic rather an elliptic problem. In fact, in 
the infinite horizon case one can obtain a closed form representation for the value function, which is not 
possible in the finite horizon case. We make use of the results of [8], and Chapter 2 of [10] (also see 
Chapter 7 of [13]) to study the properties of the approximating sequence. For example, we show that 
the approximating sequence is bounded with respect to the Holder semi-norm (see page 61 in [7] for a 
definition), which is used to argue that the limit of the approximating sequence (which is a fixed point 
of J) solves a corresponding free boundary equation. 

Somewhat similar approximation techniques to the one we employ were used to solve optimal stopping 
problems for diffusions: see e.g. [2] for perpetual optimal stopping problems with non-smooth pay-off 
functions; and [6], [5] for finite time horizon American put option pricing problems for geometric Brownian 
motion. On the other hand, [1] and [11] consider the smooth fit principle for the infinite horizon American 
put option pricing problems for one-dimensional exponential Levy processes using the fluctuation theory. 
Also see [4] for the analysis of the smooth fit principle for a multi-dimensional infinite horizon optimal 
stopping problem. 

The next two sections prepare the proof our main result, Theorems 3.1, in a sequence of lemmata and 
corollaries. In the next section, we introduce the functional operator J, that maps a given function to 
the value function of an optimal stopping problem for a geometric Brownian motion. We then analyze 
the properties of J. For example, J preserves convexity with respect to the first variable; the increase 
in the Holder semi-norm after the application of J can be controlled; J maps certain class of functions 
to the classical solutions of free boundary equations. In Section 3, we construct a sequence of functions 
that converge to the smallest fixed point of the operator J. We show that the sequence is bounded in the 
Holder norm, and satisfies certain regularity properties using results of Section 2. We eventually arrive 
at the fact that the smallest fixed point of J is equal to V. As a result the regularity properties of V 
follow. 



ON THE FINITE TIME HORIZON AMERICAN OPTIONS FOR JUMP DIFFUSIONS 3 

2. A Functional Operator and Its Properties 

Let us define an operator J through its action on a test function / : M + x M + — > M + : The operator J 
takes the function / to the value function of the following optimal stopping problem 

Jf(x,T) = sup E x \ r e - (r+A)t A-P/(5 t °,r-t)dt + e- (r+A)T (iir-5°) + l, (2.1) 

in which 

Pf(x,T-t)=[ f(xz,T -t)u(dz), x>0. (2.2) 
Jr + 

We will extend T — ► Jf(x, T) onto [0, oo] by letting 

Jf(x,oo) = lim Jf(x,T). (2.3) 

T — >oo 

Here, S = {S t °;i > 0} is the solution of 

dS t ° = + ovS?dWi, = x, (2.4) 

whose infinitesimal generator is given by 

I o o d 2 d , 

A --=r x d^ + » x Tx (2 - 5) 

In (2.1), <S[o,t] denotes the set of stopping times of S° which take values in [0, T]. Note that 

S? = sfli, (2.6) 

where 

H t = exp(((i- laAt + aWt). (2.7) 
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The next remark characterizes the optimal stopping times of (2.1) using the Snell envelope theory. 
Remark 2.1. Let us denote 

Y t := f e-( r+A ) s A • P/(S°, T - s)ds + e- {r+x ^(K - (2.8) 
Jo 

Using the strong Markov property of S°, we can determine the Snell envelope ofY as 

sup E{F T |^} = e-( A +^J/(5°,r-t)+ F e-( r + x >\Pf(S° s ,T-s)ds, te[0,T\. (2.9) 

reSt.T Jo 

Theorem D.12 in [10] implies that the stopping time 

t x := inf{t € [0, T) : & = Y t } A T = M{t £ [0, T] : J/(S°, T — t) = (K — $)+}, (2.10) 

satisfies 

J/(x,T) = j jT" e- (r+A)t A • Pf(SlT- t)dt + e - (r+A)T *(-^ " S?J + | • (2-H) 

Moreover, i/ie stopped process ( e -( r+A )(* Ar -) Jf(S? ATx , T - tAr x ) + f^ Tx e~( r+A)s A • Pf{S%T- s)ds) t > 
is a martingale. The second infimum in (2.10) is less than T because Jf(S^,0) = (K — S^) + ■ 

When f is bounded, it follows from the bounded convergence theorem that (using the results of [3] and 
arguments similar to the ones used in Corollary 7.3 in Chapter 2 of [10]) 



J/(x,oo) 



sup E x \ f T e" (r+A)t A • P/(S t °, oo)dt + e- {r+x)T (K - S°)+l . (2.12) 
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The next three lemmas on the properties of J immediately follow from the definition in (2.1). The 
first lemma states that J preserves monotonicity. 

Lemma 2.1. Let T — > f(x,T) be non- decreasing, and x — > f(x,T) be non-increasing. Then T — > 
Jf(x,T) is non- decreasing and x — > Jf(x,T) is non-increasing. 

The operator J preserves boundedness and order. 

Lemma 2.2. Let f : M + x M + — > M + be a bounded function. Then Jf is also bounded. In fact, 

0< || J/Hoc < # + —^-11/1100- (2.13) 
r + A 

Lemma 2.3. For any /i, f 2 : x M + — > M + t/mf satisfy f\(x,T) < /2(x,T), w;e /iaue i/iai Jfi(x,T) < 
Jf 2 (x,T) for all (x,T) G E+ x M + . 

As we shall see next, the operator J preserves convexity (with respect to the first variable). 

Lemma 2.4. If f : M + x]R_|_ — > M + is a convex function in its first variable, then so is Jf : M + x]R + — > M + . 

Proof. Note that J/ can be written as 

Jf(x,T)= sup e( [ T e'^+^X- Pf(xH t ,T -t)dt + e-( r+ V T (K -xH T ) + X . (2.14) 

Since f(-,T — t) is convex, so is Pf(-,T — t). As a result the integral with respect to time in (2.14) is 
also convex in x. On the other hand, note that (K — xH T ) + is also a convex function of x. Taking the 
expectation does not change the convexity with respect to x. Since the upper envelope (supremum) of 
convex functions is convex, the result follows. □ 

Remark 2.2. Since x = is an absorbing boundary for the process S°, for any f : M+ x M + — * 

J/(0, T) = sup { f e~( r+A ) s A/(0, T - s)ds + e-^ +r ^K\ 
te{o,T} [Jo J 

= max jif, £ e~( r+A)s A/(0, T - s)ds + e~ {x+r)T , T > 0. 

If we further assume f < K , then Jf(0, T) = K, T > 0. 

Lemma 2.5. Let us assume that f : x M + — ► M + zs convex in its first variable and ||/||oo < T/ien 
x — > Jf(x,t) satisfies 

\Jf(x,T)-Jf(y,T)\<\x-y\, (x,y)eR+xR + , (2.16) 

and all T > 0. 

Proof. First note that a positive convex function that is bounded above has to be non-increasing. There- 
fore / is non-increasing. As a result of Lemma 2.1, x — > Jf(x,t) is non-increasing. This function is 
convex (by Lemma 2.4), and it satisfies 

J f(x, T) > (K - x) + J/(0, T) = K. (2.17) 

Consequently, the left and right derivatives of Jf satisfy 

- 1 < D x _Jf(x, T) < D x + Jf(x, T) < 0, x > 0, T > 0. (2.18) 

Now, the result follows since the derivatives are bounded by 1 (also see Theorem 24.7 (on page 237) in 
[17]). □ 



(2.15) 
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Remark 2.3. Let Tq £ (0, oo) and denote 

F{x,T) = sup E\e-( r+x ^(K -xH T ) + \ , x G R + , T G [0, T ]. 

Then for S <T <T 

F(x, T) - F(x, S)<C-\T- S| 1/2 , 
for all x G R+ and for some C that depends only on Tq. See e.g. equation (2.4) in [14]- 



(2.19) 
(2.20) 



The next lemma, which is very crucial for our proof of the smoothness of the American option price 
for jump diffusions, shows that the increase in the Holder semi-norm that the operator J causes can be 
controlled. 



Lemma 2.6. Let us assume that for some L G (0, oo) 

\f(x,T)-f(x,S)\ <L\T-S\ l '\ (T,S) e [So, lb] x [S ,T ], 
for all x G R+, for < So < To < oo. Then 

\Jf(x,T)-Jf(x,S)\ < (aL + C) \T-S\ l '\ (T, S) G [S ,T ] x [S ,T ] 
for some a G (0, 1) whenever 

\T-S\< 

Here, C G (0, oo) is as in Remark 2.3. 



r L 
r + XXK 



Proof. Without loss of generality we will assume that T > S. Then we can write 



Jf(x,T) - Jf(x,S) < sup 



E 



sup 



+ 1 {S<t} 

X 



E A Pf ( XHh T ~ ^ + e ~ (r+X)T ( K - xH r) + 

|^ TA6 e^ r+x)t X Pf(xH t , S - t)dt + e -( r + A )(- AS )(ET - xH TAS ) + 1 
E { Jo e_(r+A)<A ( Pf ( xHt ' T — t) - Pf{xH u S - t)) dt 

J T e -(r+Vt X pf( x H t ,S - t)dt + {e^ r+ ^ T {K - xH T ) + - e^ r+ ^ s (K - xH s ) + 

< — —rL (T - S) 1 / 2 + K (e~^ s - e~^ T ) 

~ r + A v ' r+AV / 



+ sup E{(e-( r+ V T (K -xH T ) + } - E {e~( r+ ^ s (K - xH s ) + )} 



< — ^—L (T - S) 1 / 2 + A K (T - S) + e -^ s {F(H S , T - S) - F{H S , 0)) 



A 

r + \~ 
A 

r + X' 



(2.21) 
(2.22) 

(2.23) 



(2.24) 



L + c) (T - S) 1/2 + XK (T - S) 
in which F is given by (2.19). To derive the second inequality in (2.24), we use the fact that 

\Pf(xH t ,T-t)-Pf(xH t ,S-t)\ < f u{dz)\f(xzH u T -t) - f(xzH t ,S - t)\ < L \T - S^ 2 , 



(2.25) 
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which follows from the assumption in (2.21), and that 

E 1 1 {5<T} e ~( r + A )*^ Pf( x H u S - t)dt\ < Aif E e^+^dtX < ( e ~(H-A)S - e~^ T 

(2.26) 

To derive the third inequality in (2.24), we use 

e -(r+X)S _ e -(r+X)T < e -(r+X)S ^ + ^ _ ^ < ^ + A )( T _ 5) ( 2 _ 27 ) 

The last inequality in (2.24) follows from (2.20). Equation (2.22) follows from (2.24) whenever T and S 
satisfy (2.23). □ 

Let us define the continuation region and its sections by 

C J f ■= {(T, x) G (0, oo) 2 : Jf(x, T) > (K - x) + }, and C 3 J := {x G (0, oo) : Jf{T, x) > (K - x) + }, 

(2.28) 

T > 0, respectively. 

Lemma 2.7. Suppose that f : M + x ]R + — > K + is sitc/i that x — > f(x,T) is a positive convex function, 
T — > f(x,T) is non- decreasing, and ||/||oo < ^f- Then for every T > i/iere exists c J f(T) G (0, if) suc/i 
i/iai C^/ = (c*^ (T), oo). Moreover, T — > c*^ (T) is non-increasing. 

Proof. Let us first show that if x > if, then a: € dj for all T > 0. Let t £ := inf {0 < t < T : Sf < if - e}. 
Since P{0 < r e < T} > for x > if, for all T > 0, we have that 

j jT' e- (r+A)t A Pf(Sf,T- t)dt + e^ r+x ^(K - 5° ) + | > 0, (2.29) 

which implies that x G C T . On the other hand, it is clear that 

(if - x) + < Jf(x,T) < J/(x,oo), (x,T) G R+ x M + . (2.30) 

Thanks to in Lemma 2.6 of [3], there exist 1* G (0, if) such that 

Jf(x, oo) = (if - x) + , x G [0, I*]; J/(x, oo) > (if — x) + , x G (^,00). (2.31) 

Since x — > J/(x,oo) and x — > Jf(x,T), T > 0, are convex functions (from Lemma 2.2 in [3] and 
Lemma 2.4 respectively), (2.29), (2.30) and (2.31) imply that there exists a point c J f(T) G (/^,if) such 
that 

Jf{x) = (If -x)+, x G [0,c J/ (T)]; Jf(x,T) > (if - x) + , x G (c J/ (T), 00), (2.32) 

for T > 0. This proves the first statement of the Lemma. The fact that T — > c(T) is non-increasing 
follows from the fact that T — > Jf(x,T) is non-decreasing. □ 

In the following lemma we will argue that if / has certain regularity properties, then J f is the classical 
solution of a parabolic free boundary equation. 

Lemma 2.8. Let us assume that f : M + x M + — > R + is convex in its first variable, ||/||oo < K and 
T — > f(x,T) is non-increasing. Moreover, we will assume that f satisfies 

\f(x,T) - f(x,S)\ < A\T - S\ 1/2 whenever \T - S\ < B, (2.33) 

for all x G R+, where A,B are strictly positive constants that do not depend on x. Then the function 
Jf : x — > is the unique bounded solution (in the classical sense) of 

Au(x,T)-{r + X)-u(x,T) + X-(Pf)(x,T)--^u(x,T) = x>c Jf (T), (2.34) 

u(x, T) = (if — x) x<c Jf (T), (2.35) 
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in which A is as in (2.5) and c 3 * is as in Lemma 2.7. Moreover, 

AJf(x,T)-(r + X)-Jf(x,T) + \-(Pf)(x,T)-^ T Jf(x,T)<0, x<c J f{T). (2.36) 

Proof. The proof is motivated by Theorem 2.7.7 of [10]. Equation (2.35) is clearly satisfied by Jf. In 
what follows, we will first show that Jf satisfies (2.34). Let us take a point in (t,T) G C J f and consider 
a bounded rectangle R = {t\,tz) x (xi,X2) containing this point. We will let 

rA V 
(r + 

Let 8qR be the parabolic boundary of R and consider the parabolic partial differential equation 

Au(x,T)-(r + X)-u(x,T) + X-(Pf)(x,T)--^u(x,T) = inP, 

u(x,T) = Jf(x,T) on doR. 

As a result of Lemmas 2.5 and 2.6, Jf satisfies the uniform Lipschitz and Holder continuity conditions, 
which implies that Jf is continuous. On the other hand, for any (T, x) G R 

\Pf(x,T) - Pf(y,S)\ < \Pf(x,T) - Pf(x,S)\ + \Pf(x,S) - Pf(y,S)\ 



< / u(dz)(\f(xz,T)-f(xz,S)\ + \f(xz,S)-f(yz,S)\) ( 2 .39) 

<A \T-S\ 1 ' 2 + i \x-y\, 

Now, Theorem 5.2 in [8] implies that (2.38) has a unique classical solution. We will show that this unique 
solution coincides with Jf using optional sampling theorem. Let us introduce the stopping time 

r := mf{9 G [0,t - h) : (t - 9,x H e G d R} A (to - ti), (2.40) 

which is the first time 5° hits the parabolic boundary when S° starts from (xo, to)- Let us a ls° define the 
process N e := e -( r+A ) e u(x i^, t - 9) + /* e'^'+^X ■ Pf(S$, t - t)dt, 6 G [0, t - ti]. From the classical 
Ito's formula it follows that the stopped process Nq At is a bounded martingale. As a result 

u(*o,t ) = = IE" {N T } = IE L-( r+x >Jf(xH T ,t - r)) + ^V^+^A • P/(5?,t - *)dtl - (2.41) 

Clearly r < r x . Since the stopped process ( e -( r + A )( tAr -) Jf{S^ Tx , * - Mt x ) + f* AT * e"( r+A )* A • P/(S°, t - 
s)ds)t>o is a bounded martingale, another application of the optional sampling theorem yields 

IE |e-( r+A ) T Jf(x H T , t - r) + j* e^'+^X ■ P/(S t °, t - t)dt J = Jf(x , t ). (2.42) 

Combining (2.41) and (2.42), we see that (2.34) is satisfied in the classical sense since the choice of 
(xo,to) G C J f is arbitrary. 

We still need to show uniqueness among bounded functions. Fix x > c J f(T). Let a be a bounded 
function satisfying (2.34) and (2.35). Let us define M t := e~^ r+x ^u(xH t ,T -t) + jj e^ r+x ^X- Pf(S° s ,T - 
s)ds. Using the classical Ito formula it can be seen that Mtf\ Tx is a bounded martingale. Since t x is optimal 
(see (2.11)), by the optional sampling theorem, we have 



i(x,T) = M = E x {M Tx } = E ^ r+x ^u(xH Tx ,T - t x ) + £' e'^+^X ■ Pf(S° s ,T- s)ds j 
= E je-( r+A ) T (/f - xH Tx ) + + e" (r+A) *A • Pf(S° s ,T- s)ds^ = Jf(x,T). 



(2.43) 
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Next, we will prove (2.36). To this end, let x < c J f(t). Let U a closed interval centered at x such that U C 
(0,c J f(T)). Let w = {t > : xH t £ U}. Since (e"( r+A )* Jf(S§, T - 1) + jj e-^ r+x >X ■ Pf(S° s , T - s)ds) t > 
is a supermartingale we can write 



IE 



-(r+\)(TuM) 



Jf(xH TuAt ,T - Tu At))+ / e-^ u XPf(xH u , T - u)du 

Jo 



<Jf(x,T), (2.44) 



for all t > 0. Since Jf(x, t) = K — x when (T, x) € — C" 7 ^, we can apply Ito's formula to obtain that 



limE 



\ e- (r+A)u ( (A - (r + A) • --^\ J/( x fT u , t - n) + XPf(xH u , T - u)j 



du 



< 0. (2.45) 



Now, (2.36) follows thanks to dominated convergence theorem, which allows us to exchange the limit 
and the expectation. We can apply the dominated convergence theorem thanks to the fact that U is a 
compact domain. □ 

Lemma 2.9. For a given T > 0, let x — > f(x,T) be a convex and non-increasing function. Then the 
convex function x — > Jf(x,T) is of class C 1 at x = c(T), i.e., 



|j/(x,T) 



x=c(T) 



(2.46) 



Proof. The proof is similar to the proof of Lemma 7.8 on page 74 of [10], but we will provide it here for 
the sake of completeness. If we let x = c(T), then 



Jf(x + e, T) = E |^ Ta:+£ e" (r+A) * A • P/((x + e)H u T - t)dt + e -(''+ A )^ (Jf _ ( x + e)#- 

e -(r+A)i A . pf( x H t , T - t)dt + e ~ {r+Xr ^\K - xH T ) 



T~X + E , 



= E 
+ E 



Tx + E 



Tx+E 



(2.47) 



e -(r+A)t A . + £)HuT _ t) _ Pf (xH u T - t)} dt 

+ E { e ~(r+x)rx + E ^ K _ {x + £ )H Tx+s )+ — (K — xH Tx+e )+] } 

< J fix, T) + E {l {Tl+e < T} e- (r+A)T ^ [{K -(x + e)H Tx+e ) - (K - xH Tx+ J] } 
+ E [l {Tx+e=T} e-( r+x ^+< [(K - (x + e)H Tx+E ) + - (K - xH Tx+s )+] } 

< Jf(x,T) - eW {l {Tx+e < T} e- {r+X)T * +s H Tx+s } 

= Jf(x,T) - eW {e-( r+x ^H Tx+E } + eW {l {Tx+e=T} e~^ T H T ] . 

The first inequality follows since t x+£ is not optimal when S° starts at x and x — > Pf(x, T) is a decreasing 
function for any T > 0. From (2.47) it follows that 

D$Jf(x+,T)<-l, (2.48) 

since e~^ r+x ^ >t Ht is a uniformly integrable martingale and t x + £ j. 0. Convexity of Jf(t,x) (Lemma 2.4) 
implies that 

- 1 = D x LJf{x-,t) < DlJf(x+, t) < -1, (2.49) 
which yields the desired result. □ 
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3. A Sequence of Functions Approximating V 
Let us define a sequence of functions by the following iteration: 

v (x,T) = (K - x) + , v n+1 (x,T) = Jv n (x,T), n > 0, for all (x,T) G R + x R+. (3.1) 
We extend these functions onto M+ x M + by letting 

v n (x,oo)= lim v n (x,T). (3.2) 

T — >oo 

This sequence of functions is a bounded sequence as the next lemma shows. 
Corollary 3.1. For all n > 0, 

(K - x) + < v n (x, T) < (l + ^) K, (x, T)eR + x M.+. (3.3) 

Proof. The first inequality follows since it may not be optimal to stop immediately. Let us prove the 
second inequality using an induction argument: Observe that vq(x,T) = [K — x) + , (x,T) G R + x R + , 
satisfies (3.3). Let us assume that (3.3) holds for n and show that it holds for n + 1. Using (2.13), we 
get that 

K+llloc = \\JvnWco < K + 1 1 + - ) K = 1 1 + - ] . (3.4) 

r + A\ r J \ r J 

□ 

As a corollary of Lemmas 2.3 and 2.4 we can state the following corollary, whose proof can be carried 
out by induction. 

Corollary 3.2. The sequence (v n (x,T)) n >o is increasing for all (x,T) € M+ x M + . For each n, the 
function x — > v n (x, T), x > 0, is convex for all T G M + . 

Remark 3.1. Let us define, 

Voo(x,T) := supv n (x,T), (x,T) G M.+ x M+. (3.5) 

n>0 

This function is well defined as a result of (3.3) and Corollary 3.2. In fact, it is convex, because it is the 
upper envelope of convex functions, and it is bounded by the right-hand- side of (3.3). 

Corollary 3.3. For each n > andt G M+, x — > v n (x,T), is a decreasing function on [0, oo). Moreover, 
T — > v n (x,T) is non- decreasing. The same statements hold for x — > Voo(x,T), and T — > v 00 (x,T), 
respectively. 

Proof. The behaviour with respect to the first variable is a result of Corollary 3.2 and Remark 3.1 
since any positive convex function that is bounded from above is decreasing. For each n, the fact that 
T — > v n (x,T) is non-decreasing is a corollary of Lemma 2.1. On the other hand, for any T > S > 0, we 
have that v QO (x, T) = sup n v n (x, T) > sup n v n (x, S) = v QO (x, S). □ 

Next, we will sharpen the upper bound in Corollary 3.1. This improvement has some implications for 
the continuity of x — > v n (x, T), n > 1, and x — ► ^00(^5 T) at x = 0. 

Remark 3.2. The upper bound in (3.1) can be sharpened using Corollary 3.3 and Remark 2.2. Indeed, 
we have 

{K-x) + < v n {x,T) < K, for each n, and (K — x) + < v^x, T) < K, (x, T) G (0, oo) 2 . (3.6) 

It follows from this observation that for every T G M+, x — > v n (x, T), for every n, and x — ► ^00(^5 T), are 
continuous at x = since v n (0,T) = ^(OjT) = K and these functions are convex. (Note that convexity 
already guarantees continuity for x > 0.) 
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Lemma 3.1. The function Vqo is the smallest fixed point of the operator J. 
Proof. 

Voo(x,T - t) = supv n (x,T - t) 

n>l 

= sup sup E x { [ T e^ r+x)t X ■ Pv n {S?,T - t)dt + e~ {r+x)T {K - S°)+ 

= sup supE*{ r e-^ r+x ^\-Pv n {S Q u T -t)dt + e-^ r+x ^ T {K - S° T ) + \ (3.7) 

= sup E x \ [ T e-^'X- P( S npv n ){S^,T -t)dt + e~^ T (K - S^) + \ 
T€5 ,r Uo n>l J 

= JVoo(x,T - t), 

in which the fourth equality follows by applying the monotone convergence theorem three times. Let 
w : x M + — ► IR + be another fixed point of the operator J. We will argue by induction that w > v^. 
For G E+ x H + , w(x,T — t) = Jw(x,T — t), which implies that w{x,T — t) = Jw(x,T — t) > 

(K — x) + = Vq(-). If we assume that w(x,T — t) > v n (x,T — t), then w(x,T — t) = Jw(x,T — t) > 
Jv n (x,T — t) = v n+ i(x,T — t). Consequently w(x,T — t) > v n (x,T — t) for all n > 0. As a result 
w(x, T -t)> sup n > v n (x, T -t) = Voo(x, T -t). □ 

Lemma 3.2. The sequence {v n {-,-)} n >o converges uniformly to Voo. In fact, the rate of convergence is 
exponential: 

v n (x,T) < Voo(x,T) < v n (x, 

T ) + (lT7) K > 0,T)gR + xR + . (3.8) 

Proof. The first inequality follows from the definition of Voq. The second inequality can be proved by 
induction. The inequality holds when we set n = by Remark 3.2. Assume that the inequality holds for 
n > 0. Then 

v O0 (x,T)= sup E x ( [ T e-^+^X-PvooiS^T-^dt + e-^+^iK - S°) + \ 
t€S ,t I Jo J 

sup w{ I" e-t'+^X- Pv n (S?,T -t)dt + e-( r+x)T (K - + [°° dt e~( A+r )*A f-A- ^ K 

g5 ,t Uo J Jo \X + rJ 



< 



n+1 



= Vn+l (x,T)+^—) K. 

(3.9) 
□ 

Remark 3.3. Note that, for a fixed Tq > 0, 

v n (x,T) < Voo (x,T) < Vn (x,T) + (l- e -( r+A ) To )" (j^A K, x € M+, T € (0, Tq). (3.10) 
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This can be derived using an induction argument similar to the one used in the proof of Lemma 3.2. We 
simply replace (3.9) by 



T)< sup E x \ f V (r+A) 'A- Pv n {S?,T -t)dt + e- {r+x)T (K - 

tgs 0;T [Jo J 



(3.11) 

Observe that one can replace K in (3.10) by \\vqo — i>o||oo- Note that the convergence rate in (3.10) is 
fast. This will lead to a numerical scheme, whose error versus accuracy characteristics can be controlled, 
for pricing American options. 

Remark 3.4. Let To £ (0, oo). It can be shown using similar arguments to the ones used in the proof of 
Lemma 2.6 that 

\ Vl (x,T) - Vl (x,S)\ < XK \T - S\ + C\T - S\ 1/2 , T,Se(0,T ], (3.12) 

for all x £ R+, in which C € (0, oo) is as in Remark 2.3. In fact 

\vx(x,T) - Vl (x,S)\ <L\T-S\ 1/2 , (3.13) 

for all x £ R+ and for some L that depends only on To. 

The next lemma shows that the functions v n , n > 0, and -Uqo are locally Holder continuous with respect 
to the time variable. 

Lemma 3.3. Let To € (0, oo) and L € (0, oo) be as in Remark 3.4 and C € (0, oo) be as in Remark 2.3. 
Then for T,S € (0, To), we have that 

\v n (x,T)-v n (x,S)\<(L + -^—]\T-S\ 1 / 2 whenever \T - S\ <( ^— -^-) , (3.14) 

\ 1 — a J \r + A XK J 

for all x € M+ and for all n > 1. Here, a € (0, 1) is as in Lemma 2.6. Moreover, 

\v 0Q {x,T)-v 00 {x,S)\<(l + -^—]\T-S\ 1 / 2 whenever \T - S\ <( ^— -^-] , (3.15) 

\ 1 — a J \r + A XK J 

for all x € R+ . 

Proof. The proof of (3.14) will be carried out using an induction argument. Observe from Remark 3.4 
that (3.14) holds for n = 1. Let us assume that (3.14) holds for n and show that it holds for n + 1. Using 
Lemma 2.6, we have that 



L + j^j+C^j \T-S\ 1 ' 2 , 



v n+1 (x,T)-v n+1 (x,S)\ < [a [L + - )+C )\T-S\ 1 / 2 , (3.16) 



for \T — S\ < ^;q_x L+C \k ) ■ ^ * s c l ear that the right-hand-side of (3.16) is less than that of (3.14), 
and 

r i+o/a-.) 2 r L 



r + X XK ~ r + X XK' 
from which the first statement of the lemma follows. Now let us prove (3.15). To this end observe that 

\voo(x,T) - Voo(x,S)\ < \voo(x,T) - v n (x,T)\ + \v n (x,T) - v n (x,S)\ + \voo(x,S) - v n (x,S)\ 
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for any n > 1, which follows from (3.14) and Lemma 3.2. The result follows since n on the right-hand-side 
of (3.18) is arbitrary. □ 

Lemma 3.4. For n > 0, \v n (x,T) - v n (y,T)\ <\x — y\, and \voo{x,T) - v 00 (y, T)\ < \x — y\, (x,y) € 
M+ x R + , for all T > 0. 

Proof. It follows from Remark 3.2 that ||f n ||oo < K, for all n > 0, and ||uoo||oo < K. Moreover, for each 
n>0, v n (-,T) is convex (for all T £ R + ) as a result of Corollary 3.2. On the other hand, it was pointed 
in Remark 3.1 that v 00 (-,T) is convex for all T € K+. Since 

v n+ i(x,T) = Jv n (x,T) and Voo(x,T) = Jvoo(x,T), (3.19) 

the statement of the lemma follows from Lemma 2.5. □ 

Lemma 3.5. For all T > and n > 0, +1 = (c Un+1 (T) , oc) /or some c^+^T) € (0,K) and C£°° = 
(c v °° (T), oo) /or some c 1100 E (0, -fT). T/ie function u n +i is £/ie unique bounded solution (in the classical 
sense) of 

Av n+1 (x,T) - (r + A) • Un+1 (x, T) + A • (Pt/ n ) (x, T) - ^ n+1 (x, T) = 0, x > c^ 1 (T) , 
*%+i(x,T) = (if - x), x < c^+HT), 
and it satisfies 

R 

= -1, T>0. (3.21) 



^v n+1 (x,T) 



x=c" n + 1 (T) 

Moreover, Voo is the unique bounded solution (in the classical sense) of 

Av^x, T) - (r + A) • tUx, T) + A • (Pv oc )(x, T) - JUoofo T) = x > c v °° (T), 



Uoo (x,T) = (K-x) x<c v °°(T), 
and it satisfies 

■^-Voo(x,T) 

ox 

On the other hand, 



OT y ' ' v n (3.22) 



-1, T > 0. (3.23) 



AvaofaT) - (r + \)-v 00 (x,T) + X-( y Pv 00 )(x,T) - ^- Voo (x,T) < x < c u °°(T). (3.24) 

oT 

Proof. The fact that 0+ 1 = (c^ +1 , oo) and C u °° = (c v °° , oo) for some c Vn+1 G (0, K) and c"°° G (0, K) 
follows from Lemma 2.7 since the assumptions in that lemma hold thanks to Corollaries 3.2, 3.3; Re- 
marks 3.1 and 3.2; and Lemma 3.1. 

The partial differential equations (3.20), (3.22); and the inequality in (3.24) are satisfied as a corollary 
of Lemma 2.8; Corollaries 3.2 and 3.3, Remarks 3.1, 3.2; Lemmas 3.1, 3.3. 

Observe that since v n is convex (Corollary 3.2) and non-increasing (Corollary 3.3) with respect to its 
first variable, t> n +i (= Jv n ) satisfies the smooth fit condition in (3.21) as a result of Lemma 2.9. The 
smooth fit condition in (3.23) holds for Voo as a result of Lemma 2.9 since t>oo (= Jvoo) (Lemma 3.1) and 
x — > Voo(x,T) is non-increasing and convex. □ 

The next lemma will be used to verify the fact that V = Voq. The classical Ito's rule can not be applied 
to the process t —>■ v oa (St,T — t) since the function may fail to be C 2,1 at T — ► c Vo °(T). As a result, 
the semi-martingale decomposition of the process t — > v 00 (St,T — t) may contain an extra term term due 
to the local time of the process S at the free boundary. 
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Lemma 3.6. Let X = {Xt;t > 0} be a semi-martingale and b : R+ — ► R be a continuous function of 
bounded variation. Let F:lx R + — > R be a continuous function that is C 2,1 on C and D (it may not 
be necessarily C 1 ' 2 across the boundary curve b), in which 

C 4 {(x, t) € R x R + : x < b(t)}, D = {(x, t) £ R x R + : x > 6(t)}. 

T/iat is, i/iere exii too functions F l ,F 2 : R x M + — > R, that C 2,1 onKx R +) and F(x, t) = F 1 (x, t) when 
(x,t) £ C and F(x,t) = F 2 (x,t) when (x,t) £ D. Moreover, F 1 (b(t),t) = F 2 (b(t),t). Then the following 
generalization of ltd 's formula holds: 

r* i 

F{X t ,t) = F(X ,0) + / -[F t (X s _+,s) + F t (X s _-,s)]ds 
Jo 1 

+ \f [F x (X s _+,s) + F x {X s _-,s)]dX s + \ [ l {Xs ^ b{s)} F xx (X s _,s)d(X,X) c s 
z Jo z Jo , 

( i >) (3-25) 

+ Yl \F(Xs,s) - F(X s -,s) - -AX s {F x {X s „-,s) + F x (X s -+,s)}\ 

+ \ j [F X (X S .+, s) - F X (X S _-, s)] l {Xs 
z Jo 



--b(s)} dL t^ 



where L t is the local time of the semi-martingale X t — b(t) at zero (see the definition on page 216 in [16]). 

Lemma 3.6 was stated in Theorem 2.1 of [12] for continuous semimartingales. The generalization for 
the case when the underlying process is not necessarily continuous is intuitively clear and just technical, 
but we will prove it in the Appendix for the sake of completeness. We are now ready to state the main 
results. 

Theorem 3.1. The value function V is the unique bounded solution (in the classical sense) of the integro- 
partial differential equation in (3.22). Given (x,T) £ R + x R + belongs to the optimal continuation region 
if x > c v (T). Moreover, it satisfies the smooth fit condition at the optimal stopping boundary, i.e., 



9 V(x,T) 



dx 



= -l,T> 0. 

x=c v ^> (T) 



Proof. The proof is a corollary of the optional sampling theorem and the generalized Ito's formula given 
above. Let T £ (0, oo) and define 

M t = e- rt v 00 (S t ,T -t), and r x := T A inf{t £ [0, T] : S t < c v °°(T - 1)}. (3.26) 

It follows from (3.22) and the classical Ito's lemma that {iW*^ }o<t<T is a bounded P-martingale. Using 
the optional sampling theorem, one obtains 

Uoo(x,r) = M = E* {M f;c } = E x [e- r ^ Voo {Sr x ,T - t x )} = ^ [e***(K - 5 ? J+} < V(x,T). (3.27) 

In the rest of the proof we will show that Voo(x, T) > V(x,T). Since Voo satisfies the smooth fit principle 
across the free boundary, when we apply the generalized Ito's formula to v 00 (S'i,T — t), the local time 
term drops. Thanks to (3.22) and (3.24), Vcx>(St, T — t) is a positive P-supermartingale. Again, using the 
optional sampling theorem, for any r £ <So,t 

Voo (x,T) = Mo > E x {M r } = E x {e- rr w 00 (5 r ,T - r)} > E x {e~ ri \K - S T ) + ] . (3.28) 

As a result Voo(x,T) > V(x,T). □ 
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Remark 3.5. We have that 

C v ™ = {xe (0, oo) : Voo > (K - x) + } = (c v °° (T), oo). (3.29) 
On the other hand, = K — x for x < c v °° . Since V = Voo, by Theorem 3.1, it follows that 

C% = {x G (0, oo) : V > (K - x) + } = (c v °° (T), oo). (3.30) 



A. Appendix 

Proof of Lemma 3.6. As in [12] we will define Z\ = X t A b(t), = X t V 6(i), and observe that 



F(X t ,t) = F\Zl t) + F\Z l t , t) - F(b(t),t). (A.l) 

On the other hand, applying the Meyer-Ito formula (see Theorem 70 in [16]) to the semi-martingale 
Xt — b(t), we obtain 



\X t - b(t)\ = \X Q - 6(0)| + / sign(X s _ - b(s))d(X s - 6(a)) 

Jo 

+ 2 E [ l {x^>b(s)}( X s ~ &(*))" + l {Xa <b(s)}(X s - b(s))+) + L 



0<s<t 



Since Z\ = \{X t + b(t) - \X t - b(t)\) and Z? = \ [X t + b(t) + \X t - b{t)\), using (A.2), we get 

dZ\ = i {(1 - sign(AV - b{t)))dX t + (1 + sign(AV - b(t)))db(t) - dL b t } 
- [l{x t _>&(t)}(*t - KO)" + l{x t < bW }(A t - 6(i))+] , 

dZ 2 t = i {(1 + sign(X t _ - 6(t)))cUr t + (1 + sign(X t _ - b(t)))db(t) - dL b t } 

+ [hx t ->b(t)}( x t - + i{x t < 6 (*)}(^* - & (*)) + ] • 



(A.2) 



(A.3) 



(A.4) 



It follows from the dynamics of Z l , i G {1,2} that 

d(Z\Z% = (l {Xt _ <m + \l {Xt _ =bm y(X,X) c t = l {Xt _ <b{t)} d(X,X)l (A.5) 

where the second equality follows from the occupation density formula, see e.g. Corollary 1 on page 219 
of [16]. Applying the classical Ito's formula to F 1 (Zj,t) and F 2 (Z 2 ,t) and using the dynamics of Z l and 
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Z 2 , we get 

F\Zlt) = F\zlQ>)+ [ F?(Zl_,8)ds+ [ F l x (Z l s _,s)dZl + \ f F l xx {s, Z\_)d (Z 1 , Z x )\ 

Jo Jo z Jo 

+ Y, [F\ZI, S ) - F\z]_,s) - AZlFi{Zl_,a)] 

Os<t 

= F 1 (Z o 1 ,0)+ f F}{Zl_,s)ds+\ [\l- S ign(X s ^-b(s)))F^Z 1 s _,s)dX s 
Jo z Jo 

+ i J (1 + sign(X s _ - b(s)))F^(Zl_,s)db(s) (A.6) 

- £ [i{x s _> b(S )}(^-Ks))- + i { x s < b(S )}(^-^)) + ]i ? i(^ 1 -^) 

0<s<t 

- \ J' Fl{Z\_, s)dL\ + i jf l {Xs _ <b{s)} F l xx (Z l s _,s)d (X c , X c ) s 

[F\zls)-F(Z l s _,s)-/\Z l s F l x {Z]_,s)\], and 

0<s<t 



F 2 (Z 2 ,t) = F 2 (Z 2 ,0) + [ F 2 {Z 2 _,s)ds+ f F 2 {Z 2 s ^s)dZ 2 s + \ [ F 2 x (s, Z 2 _)d (Z 2 , Z 2 ^ 

Jo Jo z Jo 

+ £ [F 2 (Zl s) - F 2 (Z 2 _,s) - AZ 2 F^(Z 2 _,s)] 



Os<t 



= F 2 (Z 2 ,0)+ I F 2 (Z 2 ^s)ds + \ [\l + sign(X^-b(s)))F 2 (Z 2 _, S )dX s 
Jo z Jo 

+ (l- S ign(X s _-b(s)))F 2 (Z 2 _,s)db(s) (A.7) 
+ £ [l { x s _>^)}(^-^))~ + l { x s < b(s )}(A s -6( S )) + ]F ; 2 (Z s 2 „, s ) 

0<s<t 

- i jf ' F 2 (Z 2 _, s)dL\ + I jf l {Xs _ <b{s)} F 2 x (Z 2 _, s)d <A C , X c ) s 
+ £ [F 2 (Zls) - F(Zl,s) - AZ]F 2 (Z 2 _,s)]] . 



0<s<t 



By splitting each term to their respective values on the sets < b(s)}, {X s _ = b(s)} and {A s „ > 

b(s)}, it can be seen that the following four equations are satisfied: 



F\ZI,Q) + F 2 (Z 2 , 0) = F(X , 0) + F(b(0), 0), (A.8) 



f Fl(Z.-, 8 )da+ I F 2 {Z 2 s __,s)ds = \ [ F t (X s „+,s)+F t (X s „-,s)ds 
Jo Jo z Jo 



+ 







F t (b(s)+, s)l {Xs _ <b(s)} + - (F t (b(s)-, s) + F t (b(s)+, s)) l { x s „=b(s)} + F t (b(s)-, s)l {Xs _ >b(s)} 



ds, 
(A.9) 
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\ J (1 - sign(X s _ - b{s)))F l x {Z l s _,s)dX s + I J (1 + sign(X s _ - b(s)))F 2 (Z 2 _, s)dX s 

= \ ! [F X (X S _+, s) + F X (X S _-, s)} dX s , 
1 Jo 

i y"(l + sign(X s _ - 6( S )))Fi(Z s 1 _, S )d6( S ) + i yV - sign(A s _ - &( S )))F X 2 (Z S 2 _, S )d6(s) = 

/"* r 1 

F x (b(s)+,s)l {Xa _ <b ( s )} + g [F x (6(s)+,s) + F x (&(s)-, s)j l{x 3 _=ft(s)} + F x (b(s)-, s)l {Xs _ >b (s)} 



(A.10) 



db(s). 
(A.ll) 



On the other hand, (3.15) of [12] still holds: 

r* r i 

F(&(t),t) = F(6(0),0) + y ^(6(s)+,s)l {Xs _ <6(s)} + - [F t (b(s)-,s) + F t (b(s)+, S )l {Xa _ =b{s) 

+ F t (b(s)-,s)l {Xs _ >b(s)} ds 

ft r 1 

y F s (6(a)+, s)l {Xs _ <fe(s)} + - [F x {b(s)-, s) + s)l {Xa _ =b(s)} ] 



(A.12) 



+ F x (6(s)-,s)l {Xa _ >il{s)} 



d6(s), 



whose proof is carried out by using the uniqueness of finite measures on p-systems. 
Let us analyze the jump terms in (A. 6) and (A. 7). We will denote 

A:=- [1 {X .->«,)}(X. ~ &(«))" + hx 3 <Ks)}(Xs ~ Hs)) + ] F l x {Z l s _,s) 
+ [F\zl,s) - F(Zl_,s) - AZlF^Z^s)}] , 

B ■= [l {Xs _ >b{s)} (Xs ~ b(s))- + l {Xa < b{s)} (X s - b(s)) + ] F 2 (Z 2 _,s) 
+ £ [F 2 {Z 2 , s) - F(Z 2 _,s) - AZ 2 F 2 (Zl,s)]] . 



(A.13) 



(A.14) 



0<s<t 

Depending on the whereabouts of A s _ and X s with respect to the boundary curve b, A and B take four 
different values: 

(1) X s ^ > b{s) and Xt > b(t). In this case 

A = 0, B = F 2 (X S , s) - F 2 (X s „,s) - AX s F 2 (X s .,s), 

and 

A + B = F(X s ,s) - F(X S _, s) - AX S F X (X S _+, s). 

(2) X s _ > b(s) and X s < b(s). In this case 

A = _( 6 ( s ) _ + F 1 ^) - F 1 (6(s), s) - (X s - b(s))F^b(s),s) 

= F 1 (X s ,s)-F 1 (b(s),s), 

B = (b(s) - X s )F 2 (b(s), s) + F 2 (b(s), s) - F 2 (X s _,s) - (b(s) - X S .)F 2 (X S _, s) 
= F 2 (b(s), s) - F 2 (X s .,s) - AX s F 2 (X s .,s), and 

A + B = F(X S , s) - F(X s _,s) - AX S F X (X S +, s). 



(A.15) 
(A.16) 

(A.17) 

(A.18) 
(A.19) 
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(3) X s _ < 6(s) and X s > b(s). We have that 
A = -(X s -b(8))F*(X a -,8) + F 1 (b(s),s) -F\X.-,8) - (b(s)-X s _)F'(X s _,s) 

= F 1 (b(s) ,s)-F 1 (X s _ , a) - AXsF 1 (X s _ , s) , 

B = (X S - b(s))F 2 (b(s), s) + F 2 (X S , s) - F 2 (b(s), s) - (X s - b(s))F 2 (b(s), s) 
= F 2 (X s ,s)-F 2 (b(s),s). 
As a result 

A + B = F(X S , s) - F(X s _,s) - AX S F X (X S _-, s). 

(4) A s _ < 6(a) and X s < b(s). Clearly, 

A = F 1 {X s ,s) - F 1 (X s _,s) - AX s F^(X s _,s) and B = 0. 

As a result 

A + B = F(X s ,s) - F(X a -,s) - AX s F x (X s .,s). 

Now combining (A.l), (A.5), (A.6), (A.7), (A.8), (A.9), (A.10), (A.ll), (A.12), (A.16), (A.19), (A.22), 
(A. 24), we obtain 

1 /•* 

F(X t ,t) =F(X o ,0) + - / [F t {X s -+,s) + F t (X s --,s)]ds 

1 Jo 

+ - [ [F x (X s _+,s) + F x (X s „-,s)]dX s + - [ l {Xs _< b{s)} F xx (s,X s _)d(X,X) c s _ 
1 Jo 1 Jo 

F(X s ,s) - F(X s _,s) - AX s F x {X s _-,s)l {Xs _< b{s)} - AX s F x (s,X s _+)l {Xs _ >b{s)} 



(A.20) 
(A.21) 

(A.22) 

(A.23) 
(A.24) 



0<s<t 



+ \J^ [F 2 (Zl,s)-F 1 (Z 1 s _,s)]dL b t . 



The last term on the right-hand-side of (A. 25) can be written as 



\[ [F 2 AZl,s)-F\zl,s)] <U\ = \f Q 



F x {X s _+,s)-F x (X s _-,s) 



{x a _=b(s)}dL t , 



(A.25) 
(A.26) 



using Theorem 69 of [16]. On the other hand, the jump term in (A.25) can be written as 



0<s<t 



F{X s ,s) - F{X s _,s) - AX S F X (X S _-, s)l {Xs _< b{s)} - AX s F x (s, A s _+)l {Xs „ >fe(s)} 



E 



0<s<t 

This completes the proof. 



F(X S , s) - F(X s _,s) - -AX S [F X (X S .-, s) + F X (X S _+, s)} 



(A.27) 



□ 
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